Evidence for Line Nodes in the Superconducting Gap of Ki_ x Na x Fe2As 2 (x = 0,0.1) 
Single Crystals from Low-Temperature Specific-Heat Measurements 
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From the measurement and analysis ol the specific heat of high-quality Ki_ ;E Na :E Fe2As2 single 
crystals we establish the presence of large T 2 contributions with coefficients q sc ~ 30 m J /mol K 3 at 
low-T for both x = and 0.1. This evidences line-nodes in the superconducting gap and suggests 
d-wave superconductivity on all Fermi-surface sheets with an average gap amplitude of Ao in the 
range of 0.5 - 0.8 meV. The derived Ao and observed T c agree well with the values calculated within 
the Eliashberg-theory, adopting a spin-fluctuation mediated pairing in the intermediate coupling 
regime. 

PACS numbers: 74.25.Bt, 74.25. Dw, 74.25.Jb, 65.40.Ba 



In spite of the substantial experimental and theoreti- 
cal research efforts to elucidate the symmetry and mag- 
nitude of the superconducting order parameter for the 
Fe pnictides P, Q , important questions concerning the 
doping evolution of the superconducting gap remain un- 
solved [3, In the stoichiometric parent compounds 
nesting usually occurs between electron (el) and hole (h) 
Fermi-surface sheets (FSS) which is responsible for the 
presence of long-range spin density wave (SDW) order. 
Superconductivity (SC) emerges when the SDW order is 
suppressed by doping or external pressure 0. An s± gap 
symmetry (nodclcss gap function with opposite signs of 
the order parameters for el and h pockets) is believed to 
be realized in under- and optimally doped compounds, 
since the antiferromagnetic spin fluctuations (SF) on the 
vector Q = (tt, tt) connecting the el and h pockets remain 
strong in the vicinity of the SDW phase. The situation 
in the ovcrdoped regime is not so clear. With further 
doping, el (h) bands disappear. Therefore, the paradigm 
of the SF glue at the vector Q = (tt, tt) docs not work. 
However, SF have been found at some incommensurate 
propagation vectors This has led to several propos- 
als for the order parameters in heavily doped compounds: 
extended s, d, s + id wave 0, Thus, even from the- 

oretical perspective the situation is still controversial. 

One of the most interesting families from this point of 
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view is Bai_ x K :r Fe2As2. Superconductivity in this com- 
pound occurs at x sa 0.15 Q. T c increases with K doping 
up to T c ps 38 K for x=QA (optimally doped regime). 
In this region experiments such as angle-resolved pho- 
toemission spectroscopy (ARPES) [lj| and thermal con- 
ductivity [ll| show the absence of nodes in the gap, con- 
firming an s-wave character of the order parameter. The 
ovcrdoped regime is not well studied, so far. The mea- 
surements on polycrystallinc samples show that T c mono- 
tonically decreases with doping, reaching T c — 3.5 K for 
KFe 2 As 2 (K122) 0. Whereas the el pockets are com- 
pletely gone in K122, it was theoretically proposed that 
a change of the order-parameter symmetry to d, s+ wave 
with accidental nodes, or s + id should occur |3j, [6H8J. 

Until now, there is no consensus in the interpreta- 
tion of available experimental data in favor of one of 
the proposed order parameters. Indeed, measurements 
of thermal conductivity [l2 , 13 1 and of the London pen- 
etration depth [hi ] have been interpreted in terms of 
d-wave SC with line nodes on each FSS. In contrast, 
recent ultrahigh-resolution laser ARPES data for K122 
have been interpreted in terms of a rather specific par- 
tial nodal s-wavc SC having an unusual gap with "octet 
line nodes" on the middle FSS, an almost-zero gap on 
the outer FSS, and a nodeless gap on the inner FSS [15|. 
The specific-heat (SH) data for K122 above 0.4 K, only, 
were interpreted within a weak-coupling BCS-like mul- 
tiple SC gap scenario with line nodes Ta, 12\, adopting 
multiband- d- or s- wave nodal SC. But at the same time 
an extremely strong coupling and/or correlated regimes 
with heavy quasiparticles were suggested for the normal 
state. Moreover the authors of a subsequent SH study in 
magnetic fields 18| doubted such an interpretation and 
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denied the presence of a smaller second gap based on the 
observation of a non-intrinsic magnetic transition at low- 
T, probably due to unknown magnetic impurities. Natu- 
rally, under such circumstances a T 2 behavior, generic for 
line nodes in the SC gap , could not be observed. 

Hence, further studies using cleaner samples are neces- 
sary to clarify the symmetry of the SC order parameter 
and the magnitude of the coupling strength. These two 
issues are the main points of the present letter. 

Ki_ T Na x Fc2As2 single crystals with £=0 (K122) and 
£=0.1 (K(Na)122) were grown using KAs as a flux. 
SH data in various fields were obtained by a relaxation 
technique in a Physical Properties Measurement System 
(PPMS, Quantum Design). The magnetic ac susceptibil- 
ity as a function of temperature was also measured by use 
of the PPMS. The T-dependence of the electrical resistiv- 
ity was measured using a standard 4-probe DC technique. 
Two single crystals, K122 and K(Na)122, were selected 
for SH measurements below 0.4 K using the heat-pulse 
technique within a dilution refrigerator. 

The T-depcndence of the volume ac susceptibilities (x' 
and x") of K122 and K(Na)122 arc shown in Fig.Q] The 
sharp transition with ~ 100 % SC volume fraction con- 
firms the bulk nature of superconductivity and the high 
quality of our crystals. The T-dependence of the in-plane 
electrical resistivity evidences a drop to zero at 3.6 K for 
K122 and at 2.9 K for K(Na)122 in agreement with the 
diamagnetic onset seen in the xac data. 

Zero-field specific-heat, C, of K122 and K(Na)122 be- 
tween 0.1 and 10 K is shown in Fig. [2] A clear sharp SC 
anomaly is observed at T c w 2.75 K for K(Na)122 and at 
T c w 3.5 K for K122 in line with the resistivity and ac sus- 
ceptibility measurements. The jump height of C at T c is 
found to be AC/T C sa 40 mJ/mol K 2 for the K(Na)122. 
This value is a bit reduced from ss 46 mJ/mol K 2 for 
K122. A zoom into the low-T range is shown in the in- 
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FIG. 2: (Color online) Specific heat C of KFe2As2 and 
K .9Nao.iFe 2 As 2 plotted as C/T vs. T. Dashed lines: Fit 
in the normal state (see text and Eq. (SI) in Ref . I22I Inset: 
Zoom into the low-T region, solid lines: Fits using Eq. Q. 



set of Fig. [2] We observe a pronounced T 2 -behavior in 
the SH for both crystals which provides a direct evidence 
of line nodes in the SC gap of Ki_ 2 ,Na a; Fe2As2, at least 
for x < 0.1 (see below). Notice that no kinks and/or 
small hump anomalies, predicted for the low-T region by 
weakly coupled two- band approaches [l6, 17 1, have been 
observed in our measurements. 

The T-dependence of the upper critical field, H C 2, of 
K(Na)122 obtained from SH measurements in fields B || c 
and B 1 c [22| is summarized in Fig. [3J From the fit 
of our data use of the single-band Wcrthamcr-Hclfand- 
Hohcnbcrg (WHH) theory Q we obtain n H°^{B j| 
c)«lT and ^ a H°^(B lc)« 4.5 T. The H c2 anisotropy 



of K(Na)122 is found to be T = Hf 2 Xc /H^ ] w 5 near T c . 
This value is in a good agreement with the one obtained 
recently for K122 j24[. The inset of Fig. [3] presents the 
scaled upper critical field, ^qH^/T 2 , as a function of the 
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FIG. 1: (Color online) Temperature dependence of the com- 
plex volume ac susceptibility Xac = X + *x" measured with 
an amplitude of 5 Oe and a frequency of v= 1 kHz, Inset: The 
normalized in-plane electrical resistivity measured in zero- 
field. 
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FIG. 3: (Color online) Temperature dependence of the upper 
critical field, H C 2, Ko.gNao.iFe2As2. Solid lines: Theoretical 
curves based on the WHH model. The dashed lines depict 
the slopes of H C 2(T) at T c . Inset: The scaled upper critical 
fields for KFe 2 As 2 H and Ko.9Nao.iFe 2 As 2 vs. T/T c . 
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reduced temperature, T/T c , for K122 [24j and K(Na)122. 
The observed scaling is generally expected in the clean 
limit HI, where H c2 {0) oc T 2 . 

Now, we will turn to the low-T part of the SH as one of 
the central points of our letter. The T-dependence of the 
SH in the SC state allows to determine whether the gap 
function possesses nodes. For instance, the exponential 
vanishing of C as in conventional s-wave SC is caused by 
the finite gap in the quasiparticle spectrum. In contrast, 
in the case of gap nodes quasiparticles are generated to 
the largest extent in the vicinity of the gap nodes. In the 
case of line nodes the quasiparticle excitation spectrum 

28f . where k± 



takes the form Eh = h 



7 k i -\- v a k,, 27 



and fen are wavevectors perpendicular and parallel to the 
FSS, respectively, vf is the renormalized in-plane Fermi 
velocity at the position of the node, while «a ~ dA/hdk 
is the slope of the gap at the node associated with the 
dispersion of the quasiparticles along the FS. This leads 
to a density of states (DOS) linear in-energy. For the 
two-dimensional (2D) case it reads: 

where i denotes the sum over all nodes. Here, we would 
like to note that both materials considered are highly 
anisotropic, justifying the 2D approximation. The ratio 
of H C 2 (in-plane vs. out-off plane) according to Fig. [3] 
suggests a large mass ratio T 2 w 25. This value agrees 
with DFT calculations for the ratio of the correspond- 



ing squared plasma frequencies [24 



Eq. 

oc a s 



quadratic power-law dependence C. 

( ^7ei(mJ/mol K 2 



results in a 
T 2 , with 



a sc (mJ/mol K 3 ) w 0.283- 



(2) 



Ao(meV) 

taken from Refs. 0, [28| valid for SC in a clean-limit, 
where we assumed that the energy gap at zero temper- 
ature, Ao, is equal for all bands, while j e i is a renor- 
malized normal-state Sommcrfcld coefficient, which is an 
average value over all FSS (see the Supplement for de- 



tails |22|). In this context, x is a dimensionless factor 



taking into account the anisotropy of 1/vf at the posi- 
tion of the nodes with respect to the averaged isotropic 
value. According to our preliminary ab initio calcula- 
tions for the fc-dependent DOS near the Fermi energy 
in dependence on the nodal direction and the expected 
amount of the anisotropic mass renormalization due-to 
anisotropic SF 0.8 < x < 1.2. The upper bound is also 
in agreement with the slightly larger gap value obtained 
in our Eliashbcrg-thcory based calculations (see below) 
than the value derived from Eq. ©. 

To analyze the low-T SH data we use the ansatz: 

C = lr T + aT 2 + f3 3 T 3 , (3) 

where the cubic term comes from phonons, which is 
defined from the data in the normal state (see Supple- 
ment [22j). The remaining terms stem from a magnetic 



cluster-glass contribution reported previously in Ref. [26J 
and from electronic degrees of freedom, which is the focus 
here. The best fit using Eq. (J3j> yields j r — 43 mJ/mol 
K 2 and a = 34.9 mJ/mol K 3 for K122 and j r = 
56.1 mJ/mol K 2 and a = 30.6 mJ/mol K 3 for K(Na)122. 
It has been shown in Ref. [26| that the residual " Sommer- 
feld" coefficient j r stems mainly due to a cluster-glass 
phase. However, we emphasize that the T 2 contribution 
measured by a is not an extrinsic "glassy "contribution 
but stems from the generic electronic contribution in the 
superconducting state, only. The estimated contribu- 
tion of the cluster glass amounts «sg ~ 2 mJ/mol K 3 
for K122 and a S G ~ for K(Na)122 [52, Hence, 
we argue that our observed large T 2 terms come from 
quasiparticle excitations near line nodes. The result- 
ing coefficients are a sc ~ 32 mJ/mol K 3 for K122 and 
a sc « 31 mJ/mol K 3 for K(Na)122. Such a large T 2 
contribution to the low T SH provides evidence of the 
d-wave character of the gap function. 

Now, we will check whether the value of a sc can be 
satisfactory described by the standard e?-wave mechanism 
on all FSS. Our analysis is based on Eq. ©, which was 
obtained within a BCS-likc theory. But it works very 
well also in the case of strongly coupled cuprates [29(. 
Therefore, we expect its validity for pnictides, too. In 
addition, we analyse the relevant coupling constants with 
the help of the strong-coupling Eliashberg equations in 
the clean limit. The coupling constants A s f of the el-SF 
and A p h of the electron-phonon interaction together with 
the two bosonic spectral densities a 2 (ui)F(u) determine 
the gap value and T c as well. A broad spectral density 
for intraband SF as obtained by recent inelastic neutron 
scattering measurements (EJ can be described as: 



a 2 F(uj) = 



7T 



r- 



(4) 



where T s f = 7.9 meV. For the phonons we assume a sharp 
Lorentzian spectral density with a width of 0.5 meV 
centered at w p h ~ 20 meV. Notice that in our simple 
isotropic effective single-band model considered here for 
a d-wave superconducting order parameter, the el-ph cou- 
pling (EPC) drops out from the gap equation by symme- 
try and it enters only the equation for the Z function 
which describes the bosonic mass renormalization of the 
h-like quasiparticles [22|. Thus, here the EPC results 
in a slight suppression of T c and Ao in contrast to the 
s± case, where the superconductivity becomes stronger 
taking into account an intraband EPC. The SC gap and 
T c calculated within Eliashbcrg-theory and their depen- 
dence on the SF coupling constant A s f for X p h = for 
the sake of simplicity and also for the more realistic case, 
with a weak el-ph coupling constant X p h ~ 0.2 (accord- 
ing to density functional theory based calculations Q) 
included, are shown in Fig. S2 of Ref. |22|. In order to 
reproduce T c , A^. s f ~ 0.64 = 0.8A ZjS f should be adopted 
for K122, where the former describes the strength of the 
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paring interaction and the latter stands for the mass 
renormalization: see also Ref. [22|. (In the isotropic s- 
wave case A^ iS f = A ZiSf holds.) These values for the cou- 
pling constants together with the calculated EPC con- 
stant A p h ~ 0.2 [24j yield a total coupling constant of 
A Zi tot ~ 1- As a result we get A(0) ET «0.75 meV (see 
Ref. [22| for more details), which is close to that in BCS 
theory cs ~ 0.65 meV obtained from the 2A /T C ra- 
tio for a d-wave superconductor in the 2D case: 



2A 3CS 



4tt 



4.28, 



(5) 



where je is the Euler constant 7^ = 1.781. The obtained 
A z ,tot values suggest that we are in the regime of inter- 
mediate coupling and we may apply Eq. ^ to extract 
Aq. Considering the isotropic case (x = 1) of Eq.@, 
the same gap function for all bands and using the experi- 
mental electronic linear coefficient 7 e / =52-68 mJ/mol K 2 
[2^ |. we estimate Ag Xp in the range 0.5 - 0.6 meV, which 
is close to the values obtained above theoretically. It is 
also close to the larger gap values Ai ~ 0.7 to 0.8 meV 
obtained within effective weak-coupling two-band mod- 
els for the electronically weakly connected bands II all 71 . 
However, the eight times smaller second gap A2 [la Il7| 
with a comparable partial DOS for the second effective 
band would produce too large a sc values exceeding our 
experimental value by a factor of three or more. Hence, 
such a multiband scenario can be excluded on the basis 
of our study [3(| ■ Thus, our d-wave scenario with compa- 
rable gaps with nodes in the gap function on all hole-like 
bands is in a good agreement with the experimental a sc 
values. Obviously, our low-T data exclude the relevance 
of any full-gap scenario derived from SH data at rela- 
tively hi gh t emperature below T c , only, such as done in 
Refs. HI, 
AC/T clel 
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In this context, we note that the SH jump 
0.7 - 0.9 is close to the theoretical weak 
coupling value of 0.95 for a d-wave superconductor (20j . 

Now, we compare the experimental value of a sc with 
that for the "octet line nodes "scenario proposed in Ref. 
[l5l ] . As shown in our Supplement j22| , the experimental 
value of a sc is too large to be described by a single gap 
with "octet line nodes"on the middle (£) FSS, only, as 
suggested in Ref. |l5|. Thus, we argue that the recently 
proposed interpocket s± scenario |8j with accidental line 
nodes on a particular FSS, only, is rather unlikely to 
be realized in the bulk as probed by the SH. Therefore, 
despite the multi-band topology of Ki_. r Na :E Fe2As2, ac- 
cording to our SH data it behaves more or less like a 
single-band d-wave SC with corresponding line nodes on 
all FSS as suggested by Reid et al. based on a thermal- 
conductivity study and universal scaling arguments valid 
for single-band d-wave SC [ll], EH and predicted also by 
microscopic weak-coupling theory 

In summary, high-quality Ki_ :c Na :c Fe2As2 (x=0,0.1) 
single crystals were studied by specific-heat measure- 
ments. Large T 2 contributions at low-T evidence the 



presence of line nodes in the superconducting gap. From 
the experimental value of a sc , an effective gap ampli- 
tude Ag l:p ~ 0.5 - 0.6 meV was estimated. This Ag Xp 
and T c agree well with the gap value calculated within 
Eliashberg-theory for a one-band d-wave superconduc- 
tor implying a moderately strong electron-boson coupling 
constant A Zi tot ~ 1- This suggests that all Fermi-surface 
sheets have comparable gap amplitudes. Our data pro- 
vide direct evidence for d-wave superconductivity in a 
Fe-pnictide system. A detailed comparison with the 
cupratcs might be helpful to deepen our understanding 
of the paring mechanism in both unconventional and still 
challenging superconducting families. 
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SUPPLEMENT 



We present details of the specific heat (SH) measurements in external magnetic fields and of the analysis of data in the 
normal state. We also briefly sketch some aspects of the renormalization of the normal-state Sommerfeld coefficient and details 
of our Eliashberg-theory based calculations for T c and the superconducting order parameter Ao discussed in the main text. 
For the clean limit case an estimate of the coefficient a BC which enters the T 2 term in the low-T SH expansion is provided for 
the s±-scenario with accidental nodes proposed recently, too. 



Specific-heat of Ko.gNao.iFe2As2 (K(Na)122) in magnetic fields 



The temperature dependence of the SH of a K(Na)122 single crystal in magnetic fields B || c (a) and B \\ ab (b) 
is summarized in Fig. IS 1 1 The SC jump shifts and broadens systematically to lower temperature with increasing 
field. The shift is more pronounced for B |j c than for B _L c reflecting the large upper critical field anisotropy of 
this compound. The extracted temperature dependence of the upper critical fields for both field orientations are 
shown in the main text in Fig. 3. For the measurements in fields B parallel ab a small copper block was used to 
orient the sample. The SH of the copper block was determined in a separate measurement and subtracted from the 
data. However, we observed a small discrepancy between the data measured with that copper block (but subtracting 
its contribution afterwards) and without a copper block as can be seen in Fig. IS 1 1 which originates in the larger 
experimental error when using the copper block. Therefore, for the quantitative analysis of the SH we used only data 
measured without a copper block. 



Analysis of the specific-heat in the normal state 



To extract correctly a sc from the experimental data shown in Fig. 2 in the main text, at first, the normal state 
SH for both samples has been fitted by the expression accounting for the electronic, lattice, and magnetic (due to a 



cluster glass (CG) phase) contributions as proposed in Ref. IS1 



C p =y n T + ec G T'+f3 3 T d + p 5 T b , (SI) 

where j n = 95(2) mJ/mol-K 2 for K(Na)122 and 89(3) mJ/mol-K 2 for KFc 2 As 2 (K122) are linear contributions 
to the SH (denoted as nominal Sommerfeld coefficients) which are the sum of the standard intrinsic electronic 
contributions 7 e ; related to the itinerant charge carriers (quasiparticles) (the so-called Sommerfeld coefficient) and 
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T(K) 



FIG. SI: Temperature dependence of the specific heat of Ko.9Nao.iFe2As2 
axis (a) and parallel to the ab plane (b). 



in various applied magnetic fields parallel to the c 



another linear, somewhat unusual glassy magnetic contribution ^cg related to the CG-phasc 
quadratic contribution due to the CG-phase. It was found that £cg ~ for K(Na)122 and £cg : 



IS 11 ] . and £cg is a 
a 2 mJ/mol-K 3 for 

K122 [SI). Finally, the lattice contributions arc 3 = 0.556 mJ/mol-K 4 and fi 5 = 1.16- 10~ 3 mJ/mol-K 6 for K(Na)122 



and /3 3 = 0.589 mJ/mol-K 4 and f3 5 = 1.2 • 10~ 3 mJ/mol-K 6 for K122. 



Now, we briefly comment on the renormalization of the Sommerfeld coefficient due to the electron-boson interac- 
tion. The Eliashberg-theory, which we use to calculate T c and A, considers the low-energy sector of characteristic 
bosonic energies. Within this approach the relation between the SH of the "bare" electrons 7^ and the renormalized 
quasiparticles due to electron-boson interaction j e i reads: 



7oi = TcH 1 + -Wot), 



where A z .tot = ^z.sf + ^z.ph consists of two contributions: the electron-spin-fluctuation interaction X z ,sf and the 
electron-phonon interaction A z . p h (see also the next subsection.) This renormalization gives a factor of about two for 
the enhancement of the Sommerfeld coefficient and similarly for the mass in the plasma frequency which determines the 
penetration depth and the superconducting condensate density for T — > 0. The enhancement factor of the order 2.5 — 3 



of the Sommerfeld coefficient as compared to the LDA calculations "i^i DA 



10.5 mJ/mol K 3 and 7^ - 30 mJ/mol 



K 3 is due to a high-energy renormalization of the bands by self-energy effects beyond the LDA based approach seen 
also in the ARPES data and similarly in the optical mass which enters the Drude weight, i.e., the intraband plasma 
frequency at high temperature and ener gies as c omp ared to the low-energy bosons included in our Eliashberg-theory 
based approach (for details see Refs. IS2l |S3L IS4|. |S5T) . 
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Details for the d-wave Eliashberg-theory based calculations 

We consider here, for the sake of simplicity, a two-dimensional (2D) single-band <i-wave superconductor with a 
cylindrical Fermi surface parameterized by the angle 9. In general, the momentum dependence of the electron-boson 
spectral function for a given mode can be expanded in terms of Fermi surface harmonics 

a 2 F(oj,9,9>) = Y,Yj(p)c?F(w, J, J')Yj,(6'), (S2) 

,7,7' 

where the functions Yj denote the spherical harmonics. For simplicity, we assume that a 2 F(oj,9,9') is dominated by 
contributions from J = (Yb = Y Z = 1) and J = 2 (Y2 =Y$ = \/2cos(20)) and that it is diagonal in the basis J, J'. 
This is the minimal set required to model c?-wave superconductivity. Then the electron-boson spectral function can 
be parameterized via dimensionlcss coupling constants X z and A^: 

a 2 F(v, 9, 9') = ^[A 4 , z + 2X^ cos(20) cos(2#')]^) , (S3) 

i 

where i is a mode index, Fi{v) = Bi(oj)/ [2 Bi(v)dvjv\ is a normalized electron-boson spectral function, and Bi(ui) 
is the bosonic density of states. We will consider contributions from two types - a single dispersionless optical-phonon 
(i = ph) branch and a spectrum of antiferromagnetic spin fluctuations (i = sf). 
It is worth to introduce: 



K,<t>{n - m) = Xz,4>y2 I 



2vF l (v)dv 

V 2 + (U} n ~ U! m y 



(S4) 



where lo„ and u rn arc fcrmion-Matsubara frequencies. 

The set of the imaginary-axis Eliashberg-equations for a d-wave superconductor reads: 



Af \ n Sro\ ( \ (j)(iu m ) cos 2 (29) 

P „ \ \J^in Z i^rn ) + T (lOJ m ) COS 2 (29) 



and 



o)„Z (w n ) = ^ + ^ A 8 (m - n) ( 7 == =P=|fP= ; 2 , 0/n ) » ( S6 ) 
P „ \ V<2 2 (!w m ) + 2 (iw m ) cos 2 (26>) / g 

where j3 = 1/keT is the inverse temperature and (■ ■ - )e denotes a Fermi surface average. The gap function is then given 
by A(iu> n ) = 4>(i(jj n )/Z(iuJn)- The transition temperature T c is determined from the highest temperature at which <f> 
has a non-zero solution while the gap magnitude A(T= 0) is approximated with a low-temperature (T ~ 0.2 K) value 
of A(ioj n = in//3). 

To describe the spin fluctuations we adopt the usual form for the bosonic density of states B s f(v) = v 2+y' 2 1 

with T s f = 7.9 meV. For the optical- phonons we assume a narrow Lorcntzian line-shape eB p h(v) = . _n Ph \2 T + r'i 

, ,^ ph /^, r i , centered at fl p h = 20 meV and T p h = 0.5 meV. We further assume that the phonons do not contribute 

[y+^Lph) +1 ph r y 

to the d-w&ve pairing (\tf>,ph = 0) and we adopt the so called unbalanced scenario for c?-wave supercondutivity [S7], 



<V, S / = 0.8A ZjS f . This ratio is also in good agreement with \ s c i in Ref. |S8 . 

In Fig. [S2]we plot 2A(r=0)/keT c extracted by use of this model. For A s f and A p h = 0, this ratio tends to about 
4.28 (not shown in Fig. S2). Thus, the electron-phonon coupling reduces the strong-coupling corrections in this special 
case. 



Comparison of a sc for d— and s±-wave (with accidental nodes) SC order parameters 

The experimental value of a S c discussed in the main text is too large to be described by a gap with "octet-line 
nodes" on the middle (£) FSS, only, as suggested in Ref. Isil. In fact, let us consider the ratio between a^ c for a d-wave 
gap with 16 nodes on all FSS having nearly equal gaps Aq ~ 0.6mcV and a| c for a single s-wave gap with 8 nodes on 
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a particular FSS. Using the general expression for the coefficient a sc in the T 2 term entering Eqs. (f ) and (2) in the 
main text, we have: 



ih 



Here, we assumed that dA/dk = dA/dO-^ and then N l (ep) is the DOS in the normal state corresponding to the i 
FSS with a node. According to the considerations given in the main text we have: 

asc * \ ai)kl dKJde ? " iJvi(£p) * °- 283 ^- (S8) 

The last expression is written for the case of a d-wave SC gap (dA/d9 = 2A ) measured in meV and where ( denotes 
Ricmann's zeta function and C(3) = 1.2. We used also that the bare Sommcrfcld coefficient 7 c i = T 3 B N to t{£F)- For the 
sake of simplicity, we also assumed that all FSS have similar k values, therefore, J^i x 1 N' l (e-p) = AxN l (e-p) = xN to t(0)- 
Finally, from Eq. jS8} for T < Af in ps 1.3 K we estimate: 

a s sc ~ 2^C7VC( £F )(l + Af)/9A s /a0 ~ ' 1 j 

where we used (according to our calculations and Ref. (s3 |) that ~ 5 with A^ ps A tot , x/x^ ~ 1.3, 

and dA s /dO pa 1.8<9Arf/c)fl, where >r^ and dA s /dO have been calculated at the nodal directions of the £ FSS according 
to the data given in Ref. |S9| . For the physical consequences of this estimated ratio with respect to the symmetry of 
the superconducting order parameter, see the main text. 
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